Alternative pathways of dewetting for a thin liquid two-layer film 
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Abstract 

We consider two stacked ultra-thin layers of different liquids on a solid substrate. Using long-wave 
theory, we derive coupled evolution equations for the free liquid-liquid and liquid-gas interfaces. 
Depending on the long-range van-der-Waals forces and the ratio of the layer thicknesses, the system 
follows different pathways of dewetting. The instability may be driven by varicose or zigzag modes 
and leads to film rupture either at the liquid-gas interface or at the substrate. We predict that the 
faster layer drives the evolution and may accelerate the rupture of the slower layer by orders of 
magnitude promoting thereby the rupture of rather thick films. 

PACS numbers: 68.15.+e, 81.16.Rf, 68.55.-a, 47.20.Ky 
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Instability phenomena in ultra-thin soft matter films with thicknesses below 100 nm be- 
came relevant mainly because they obstruct the fabrication of homogeneous coatings l|. 
The interest was further boosted by the possibility to control such processes and, to use 
them to manufacture functional layers on the nanometer scale . The stability of ultra- 
thin films is dominated by the effective molecular interactions between the substrate and 
the hhn surface Q. TW represent, for instance, long-range van-der-Wa* forces which 



5]. 



However, 



increase (decrease) the pressure in the film if they are attractive (repulsive) 
to determine the emerging length scale and pattern for unstable films a study of the film 
dynamics is required. Using a film thickness evolution equation obtained by long-wave ap- 
proximation .6? , the dewetting of a single layer of liquid is now reasonably well understood 
(see e.g. Ref. |7|). 

However, little is known on the behaviour of two stacked ultra-thin layers of simple or 
polymeric liquid on a solid substrate (see Fig.Q). Such a two-layer film allows for richer 
dynamics than an one-layer system because, both, the free liquid-liquid and the free liquid- 
gas interface evolve in a coupled way. The evolution is driven by the effective molecular 
interactions between all the three interfaces separating the four material layers: substrate, 
liquidi, liquid2 and ambient gas. Although experimental studies investigated different as- 
pects of dewettin g fo r two-layer films like interface instabilities or the growth of holes 
El 1Q, U, 12, 13, 14, ^| up to now no general theoretical description of the inter- 



face dynamics has been given [16]. The case of small interface deflections was investigated 
in Ref. for a thickness of the lower layer, di, much larger than that of the upper one 
(d 2 - d x ). 

The most intricate questions for the first stage of dewetting of a two-layer system are 
which interface will become unstable, where does the film rupture, and how long will it 
take. This will determine the observability of the instability and the final morphology of the 
film. Expe.hne.nts found roughening of the lio.uid-lio.uid interface Q or an instaMity of 
the liquid-gas interface (a, UM ■ Holes that evolve solely in the upper layer were also studied 



In this Rapid Communication, we derive and analyse coupled long-wave evolution equa- 
tions for the two interfaces that are valid for all interface deflections and thickness ratios. We 
show that solely by changing the thickness ratio of the layers one switches between different 
dominant instability modes. This leads to drastic changes of the pathway of dewetting from 
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rupture at the substrate to rupture at the liquid-liquid interface (see Fig. El below). Remark- 
ably, for systems composed of two layers of very different thickness, i.e. with very different 
time scales for the rupture of the individual layers, the faster layer drives the evolution and 
accelerates the growth of surface modulations of the slower layer by orders of magnitude. We 
illustrate our results for two-layer systems of polystyrene (PS) and polymethylmetacrylate 
fPMMA) with silicon (Si) or silicon oxide (SiO) as substrate like studied experimentally in 
Refs. BHQ- 

We believe, our model can be used not only for the description of two-layer experiments 
with simple or polymeric liquids, but using an appropriate free energy functional also for a 
liquid film on a substrate with a stable but soft coating like a polymer brush [18] and, includ- 
ing driving terms, for studies of the transport of liquid droplets in liquid-liquid microfluidic 
systems jl9(. 

Coupled film thickness equations: We obtain evolution equations for the film thicknesses 
hi and h 2 by simplifying the Navier-Stokes equations employing long-wave approximation 
p]. Thereby a no-slip condition at the substrate, the continuity of the velocity field and the 
balance of the stress-tensors at the liquid-liquid and liquid-gas interfaces are used. Consid- 
ering an isothermal two-layer system where both layer thicknesses are smaller than 100 nm, 
we neglect gravity and solely focus on the effective molecular interaction. For simplicity we 
only regard non-retarded long-range van-der-Waals forces resulting from dipole-dipole inter- 
actions between apolar materials. However, inclusion of other forces, like e.g. short-range 
polar forces |20j|, or of slip boundary conditions [6] (that may be necessary for polymer 
films) is straightforward as for one-layer films. Details of the derivation will be presented 
elsewhere. We obtain 

where 5F/5hi with i — 1, 2 denotes functional derivatives of the total energy of the system 

F = J[p s + p vw }dx. (2) 

It contains the densities of the surface energy p s = ^[o~i(V hi) 2 + o- 2 (V h 2 ) 2 }, and of the energy 
for the van-der-Waals interaction pvw = — A g2 i s / {I2nh\) — Ai\ s j \12-nh\) — Ai 2g /[127i(h 2 — 
hi) 2 ]. The surface tensions o~\ and a 2 belong to the liquid-liquid and liquid-gas interface, 
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respectively. A g2 i s , A 2 \ s and A\ 2g are four- and three-index Hamaker constants, with sub 
scripts s, 1, 2 and g refering to the substrate, liquidi, liquid2 and gas, respectively 
The symmetric matrix of the positive mobility factors Qik reads 

Q--L( " ? I (3) 

where \i\ and /i2 are the viscosities of liquidi and liquid2, respectively. Note, that for 
d 2 — d\ <C d\ and for small surface deflections Eqs. (HJ) simplify to those of Ref. 17]. Assuming 
two identical liquids, Eqs. Q reduce to the well known one-layer equation [6 . 

To compare to the well understood one-layer systems we non-dimensionalize Eqs. (fT|) 
using scales derived from the upper layer as an effective one-layer system. We scale x with 



A up = ^7i(d 2 - d 1 ) 2 y/Tra 2 /\A 1 2g\, hi with d 2 - d\ and t with r up = A8tt 2 n 2 a 2 (d 2 - dif / A\ 2g . 
The corresponding energy scale is \Ai 2g \/lQit z (d 2 — d\) 2 . The ratios of the mean thicknesses, 
surface tensions and viscosities are d = d 2 /di, a = cr 2 /ai and \i = [i 2 /[ii, respectively. To 
compare with the lower layer as effective one-layer system one introduces in an analogous 
way the length scale Ai ow and time scale t\ ow . 

We simulate the coupled time evolution of h\ and h 2 , Eqs. in an one-dimensional 
domain using a semi-implicit time integration scheme and periodic boundary conditions. 
Initial conditions consist of flat layers with an imposed noise of amplitude 0.001. Alternative 
pathways of dewetting that occur for different thickness ratios d are presented in Fig. El using 
a Si/PMMA/PS/air system as an example. Fig.|21(a) shows that for a relatively small d = 1.4 
the two interfaces start to evolve deflections that are in anti-phase indicating the dominance 
of a varicose mode. When the liquid-gas interface approaches the liquid-liquid interface the 
latter starts to move downwards due to dynamical effects. This pathway leads to rupture 
of the upper layer, i.e. at the liquid-gas interface. On the contrary, Fig.El(b) shows that for 
a larger d = 2.4 the growing deflections of the two interfaces are in phase indicating the 
dominance of a zigzag mode. As a consequence, here the lower layer ruptures, i.e. rupture 
occurs at the substrate. 

Note, that in both, Fig.|21(a) and (b), at the moment of rupture the respective non- 
ruptured layer is also in an advanced stage of its evolution leading to subsequent rup- 
ture. This is remarkable because their time scales as effective one-layer systems are 15 
times (Fig.|21(a)) and 35 times (Fig. 121(b)) slower than the time scales for the respective 
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fast layer. The ratio of the time scales t up /t\ ow is proportional to (d — l) 5 , i.e. for a lower 
layer ten times thicker than the upper one the rupture time of the lower layer is about five 
orders of magnitude larger than the one of the upper layer. However, a simulation for a 
Si/PMMA/PS/air system with d\ = 10 and d 2 — d\ — 1 shows that at rupture of the upper 
layer at t = 0.61r up = 3.99 x 10 _5 r low the lower layer already evolved a depression of one 
fourth of its thickness. If the lower layer is the fast one, the effect also exists but is less 
pronounced. 

In both cases, the acceleration of the rupture of the slower layer is caused by the direct 
coupling of the layers via the liquid-liquid interface. The fast evolution of the thinner layer 
deforms the interface and brings the thicker layer beyond the slow linear stage of its evolution. 
If the upper layer is the driving layer the process is in addition dynamically enforced because 
the liquid-liquid interface is 'pushed away' by the advancing liquid-gas interface. 

Linear stability analysis: Deeper understanding of the different pathways can be reached 
by studying the linear stability of the initial flat layers. We linearize Eqs. Q for small 
disturbances \i ex P (Pt) cos(foc) for i = 1,2, where k, [3 and x — (Xi^X?) are the wave 
number, growth rate and amplitudes of the disturbance, respectively. The dispersion relation 
(3(k) is obtained solving the resulting eigenvalue problem (k 2 Q ■ E(k) + f3I)\ — 0, where 
Q is the scaled mobility matrix and E is derived from the free energy Eq. (j2J) as Eij = 
dhihjPvw + 3ij&ik 2 (5"i — 1, a 2 — cr, and 5ij = 1 for i = j and zero otherwise). The 
stability threshold shown as solid line in Fig.El is determined by E for disturbances of 
infinite wavelength, i.e. k = 0. The system is linearly stable for 

det E > and E n > at k = 0. (4) 

An instability sets in if at least one of the conditions (@J) is violated. This implies, that the 
two-layer film can be unstable even if both, dh^pww and dh 2 h 2 pww are positive, i.e. if the 
effective one-layer systems related to these terms are both stable. 

Fixing the Hamaker constants, i.e. the combination of materials, and changing d one finds 
a line (trajectory) in the stability diagram Fig-EI as shown for a variety of experimentally 



studied systems. Interestingly, for van-der-Waals interactions calculated as detailed in [21 1 
one can show that such a trajectory can not cross the stability threshold, i.e. it is not 
possible to stabilize such a two-layer system by solely changing d. For instance, for the 
Si/PMMA/PS/air system the second condition in (j3J) is violated for all d and the system is 



always unstable. At d = 1, i.e. for a vanishing upper layer, the system is on the boundary 
between the one- and the two-mode regions (dashed line). For 1 < d < 2.3 the unstable 
mode is an asymmetric varicose mode. A corresponding dispersion relation (3{k) is shown 
for d = 1.4 in Fig.0] (cp. the time evolution in Fig.|21(a)). For d > 2.3, i.e. for smaller 
thicknesses of the lower layer, the unstable mode is an asymmetric zigzag mode. Fig.|3]gives 
(3{k) for d = 2.4 corresponding to the time evolution shown in Fig.|21(b). For the fastest 
mode the zigzag mode is strongly asymmetric, i.e. the deflection of the liquid-liquid interface 
dominates the linear stage of the evolution. Note, that f3(k) and the dominant mode type 
depend on a and /i, whereas the stability does not. 

Further on, the simultaneous action of the van-der-Waals forces between the three inter- 
faces allows for dispersion relations with two maxima. An experimental system showing this 
unusual form of f3(k) can be realized with a substrate that is less polarisable than both lay- 
ers. This is the case for the SiO/PMMA/PS/air system 211 ] . A dispersion relation showing 
maxima of equal height is given for d = 2.16 and a = 10 in Fig.HJ The maxima at small and 
large k correspond to strongly asymmetric zigzag and varicose modes, respectively. This 
implies that the larger (smaller) wavelength will predominantly be seen at the liquid-gas 
(liquid-liquid) interface (Fig.|5](a)). Increasing (decreasing) the ratio of the surface tensions 
strengthens the smaller (larger) wavelength. This implies that solely changing a by adding 
an otherwise passive surfactant one can switch from an evolution entirely dominated by the 
liquid-liquid interface to one dominated by the liquid-gas interface. This illustrates Fig. El 
by single snapshots from the non-linear time evolutions for different a. 

To conclude, we have derived coupled evolution equations for a thin liquid two-layer film 
driven by long-range van-der-Waals forces. The system represents the most general form 
of coupled evolution equations for two conserved order parameter fields in a relaxational 
situation and is apt to describe a broad variety of experimentally studied two-layer systems 
Refs. Linear and non-linear analysis have shown that the mobilities have no influence 
on the stability threshold, but determine the length and time scales of the dynamics. We 
have shown that for a two-layer system both interface deflection modes - zigzag and varicose 
mode - may be unstable and lead to rupture at the substrate or the liquid-liquid interface. 
Remarkably, the faster layer accelerates the evolution of the slower layer even if the latter 
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is rather thick implying that its rupture time may be shortened by orders of magnitude. 
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FIG. 1: Geometry of the two-layer system. The mean film thicknesses of the lower and upper layer 
are d\ and d,2 — d\, respectively. 

FIG. 2: Snapshots from time evolutions of a two-layer film for a Si/PMMA/PS/air system at 
dimensionless times (in units of r up ) as shown in the insets, (a) At d = 1.4 a varicose mode evolves 
leading to rupture of the upper layer at the liquid-liquid interface. The ratio of the time scales 
derived from upper and lower effective one-layer system is t up /t\ ow = 0.066. (b) At d = 2.4 a 
zigzag mode evolves and rupture of the lower layer occurs at the substrate (r up /ri ow = 34.98). The 
domain lengths are 5 times the corresponding fastest unstable wave length and fj, = a = 1. 

FIG. 3: Stability diagram for fixed scaled coupling dh 1 h 2 PVW = 8Tr 2 Ai2 g /\Ai2g\. Shown are the 
stability threshold (solid line) and the boundary between unstable one-mode and two-mode re- 
gions (dashed line). The thin lines represent the trajectories for commonly studied systems: (1) 
Si/PMMA/PS/air, (2) SiO/PMMA/PS/air, (3) SiO/PS/PDMS/air, (4) Si/PS/PDMS/air, and 
(5) Si/PDMS/PS/air. The Hamaker constants were calculated as detailed below [21 1 . 

FIG. 4: Dispersion relations (3{k) for a Si/PMMA/PS/air system at d = 1.4 (dashed line, varicose 
mode) and d = 2.4 (dot-dashed line, zigzag mode) for parameters as in Fig. [2 and 4/3 (/c) for a 
SiO/PMMA/PS/air system at d = 2.16 (solid line) for fj, = 1 and a = 10. 

FIG. 5: Single snapshots from time evolutions of a SiO/PMMA/PS/air system for d = 2.16, fx = 1 
and different a. (a) The respective evolutions of the two interfaces are dominated by modes of 
different wavelength {a = 10, t = 3.83). In (b) and (c) the evolution is dominated by the liquid-gas 
and the liquid-liquid interface, respectively (a = 5,t = 3.9 and a = 100, t = 0.52). 
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